Disorder induced phase segregation in La2/'iCaii^MnO-i manganites 
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Neutron powder diffraction experiments on La2/3Cai/^AIn03 over a broad temperature range above 
and below the metal-insulator transition have been analyzed beyond the Rietveld average approach 
by use of Reverse Monte Carlo modelling. This approach allows the calculation of atomic pair dis- 
tribution functions and spin correlation functions constrained to describe the observed Bragg and 
diffuse nuclear and magnetic scattering. The results evidence phase separation within a paramag- 
netic matrix into ferro and antiferromagnetic domains correlated to anistropic lattice distortions in 
the vicinity of the metal-insulator transition. 



PACS: 75.30.Vn,73.40.-c 

Rare earth (R) doped manganites with general formula 
Ri-xAxMnOj, where A is an alkaline earth cation, show 
an intricate phase diagram with a broad phenomenol- 
ogy that ranges from ferromagnetic (FM) to antiferro- 
magnetic (AF) interactions and from metallic or semi- 
conductor behaviour to charge ordered insulators, de- 
pending wooii the level of doping and some structural 
parameterad. However, recent theoretical developments 
start to draw a unified picture that confers a prominent 
role to the intrinsic inhomogeinities present in these sys- 
tems leading^to a strong tendency towards phase sep- 
aration (PS)cl. The particular transport and magnetic 
properties of a given composition would be determined 
by the specific equilibrium sustained between competing 
segretated domains: typically ferromagnetic (FM) metal- 
lic, antiferromagnetic (AF) orbital or charge ordered, and 
paramagetic (PM) insulating patchesH 

Large differences observed in the ground state of man- 
ganites are related to the existence of two different mech- 
anisms leading to the phases mixture: Firstly, electronic 
phase separation would dominate in the extreme dop- 
ing regimes and would render nanometer scale clusters. 
Secondly, equal electronic density phase separation in- 
duced by disorder would configure a mixed phase scenario 
around a first order metal-insulator (MI) transition in the 
intermediate hole dopping range. Here the level of disor- 
der is the chief variable defining the length scale of each 
phase. Within the latter context, the MI transition is 
treated in a percolative framework where the conduction 
through metallic or insulating phases alternatively dot 
nate the transport below and above Tmi, respectively! 
Magnetoresistive properties are also explained as the re- 
sult of altering the initial equilibrium between phases, 
the applied magnetic field reinforcing and enlarging the 
FM and metallic domainSj^ the expense of the insulating 
ones (either AF or PM)BElla. 

La2/3Cai/r^Mn03 (LCMO) appears to be a good re- 
alization of a system where disordered inducedrPS may 
take place. It shows a PM to FM first orderQ transi- 
tion around Tc= 270 K, about the same temperature 
where the electronic transport undergoes a MI transi- 



tion (see Fig. Nonhomogenous states have been in- 
ferred in LCMO from cxpcj'ixncnts supporting polaronic 
descriptions of the systerrfl'Ll but lacking a microscopic 
description of the inhomogenities. Scanning tunneling 
spectroscopy has revealed the mixed phase tendencies 
of LCMO thin filmal2l but, recently, nanometer PS in 
thin films has been associated to extrinsic effects, such 
as the existence of microstrains at the interface with 
the substrateM. Nanometer regions of charge order- 
ing have been also ohsjerved by electron microdiffraction 
techniques in LCMCO and local lattice distortions asso- 
ciated to polarons have been reported frem diffuse x-ray 
experiments in the PM phase of LCMCO. 

This is consistent with Neutron Powder Diffraction 
(NPD) experiments analized in terms of pair distribution 
functions that report on local lattice distorsions, com- 
patible with lattice polarpii formation as well as with a 
mixed phase escenariolixE£l. These groups, however, do 
not incorporate in their analysis the diffuse magnetic con- 
tribution to the scattering and, consequently, fail to give 
a microscopic magnetic counterpart associated to the ob- 
served lattice distortion. The existence of two such seg- 
regated magnetic piiases (FM and PM) in LCMO has in- 
deed been reportedtZI by inelastic magnetic neutron scat- 
tering experiments but, due to the low neutron momen- 
tum transfer range explored, these authors cannot pro- 
vide evidence of the coupling of the magnetic system to 
the lattice. In this letter we provide unambiguous exper- 
imental evidence of PS at a microscopic scale in LCMO 
and follow the temperature evolution of the interplay of 
magnetic and structural degrees of freedom from 15 K 
up to 325 K. rOur results are compatible with theoretical 
predictionscTEl. 

We have combined Neutron Powder Eliffraction tech- 
niques and Reverse Monte Carlo (RMC)E£l modelling of 
the Bragg (nuclear and magnetic) and diffuse scattering 
intensities. In RMC modelling, atomic and spin coor- 
dinates for a set of ions are arranged non-periodically 
within a simulation box and varied at random without 
any intervining potential. An optimal configuration for 
the nuclear and spin variables is then found whose total 
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(nuclear and magnetic) calculated powder structure fac- 
tor, S{Q), fits the observed one in a least-squares sense. 
Within this fitted configuration, calculation of nuclear 
pair correlations functions, g{r), and spin-spin spatial 
(/i(0).^(r)} correlation functions is possible along spec- 
ified directions and planes in the crystal. 

PolycrystallineJ^a2/3Cai/3Mri03 was prepared by the 
citrate techniquetJ and characterized by x-ray diffrac- 
tion to assess the phase purity. The Curie temperature, 
Tc, was determined using SQUID magnetometry to be 
around 270 K as seen in Fig. 0. NPD data were mea- 
sured on the SLAD diffractometer at the NFL Studsvik, 
in the temperature range 15K < T < and momen- 

tum transfer Q range 1 < Q{A^^) < 10. Careful measure- 
ments of the empty container, background and proper 
substraction of these contributions as well as those from 
multiple scattering were done. Rietveld refinements of 
the lattice constants were done for a Pnma cell at all 
measured temperatures and the derived atomic coordi- 
nates used as starting configuration for the RMC runs at 
each temperature. 

For the RMC simulations we used a box equivalent 
to 6x4x6 Pnma {b long axis) single cells for which the 
La and Ca atoms are treated from the point of view of 
neutron scattering as a weighted average atom. Classical 
spins are attached to the positions occupied by magnetic 
ions whose orientation is also randomly tried. The scat- 
tering contribution from this magnetic simulation box is 
added to the nuclear one in order to make the compar- 
ison with the measured structure factor. Only Mn ions 
carry classical spins whose form factor corresponds to a 
weighted Mn'^+/Mn^+ average. This is equivalent to as- 
sume equal electronic pdensity for the various phases as 
stated in the theorjU'El. For each temperature, all nu- 
clear and magnetic cells are initially dressed with atoms 
and spins as obtained from Rietveld refinement of the 
experimental data. Minimum interatomic distances and 
total magnetjic moment constrains are applied. The RM- 
CPOW codala also provides for the adequate polycrys- 
talline averaging. 

Figure shows a representative set of experimental 
data along with the corresponding RMC and Rietveld 
fittings. Notice that RMC not only shows a better agree- 
ment than the Rietveld refinement but also provides a 
microscopic decription of the system free from "a pri- 
oristic" assumptions, as seen in Figure |^. In order to 
stress the sensitivity of our RMC models to the local 
magnetic correlations, we have calculated the magnetic 
diffuse scattering for a model in which the magnetic mo- 
ments have been randomly reassigned to different nuclear 
sites. This model therefore has the same magnetic long- 
range order (i.e. Bragg peaks), but different short range 
order, the resulting diffuse scattering exhibit magnetic 
intensity in the wrong places as shown Fig. ^jc. 

We have computed g{r)MnMn for the fitted RMC con- 
figurations for Mn-Mn pairs both in the (a,c)-plane and 
along the long axis h (Fig. ||a ). These planes and di- 



rections are treated as volumes with adjustable dimen- 
sions in the 3-dimensional configuration space. Fig. 
shows the normalized contributions of the in plane (a, c) 
and out of plane Mn-Mn to g{r)MnMn- It can be seen 
that nearest neighbours (NN) along different directions 
exhibit distance distributions centered around the same 
values (in agreement with Rietveld averages). However, 
the atom pairs in the (a,c) plane seem to account for 
the tails of the NN distance contributions, while pairs of 
atoms along the long axis have a narrower distribution. 
It must be concluded, therefore, that not only distortions 
with respect to the averaged structure exist but these are 
basically localized within the (a,c)-plane. 

Similarly, we have calculated the magnetic moment 
pair correlation function, (/x(0)./i(r)). As expected, the 
long distance asymptotic limit of these correlations fol- 
low the macroscopically measured magnetization. That 
is, above Tc the system is a paramagnet ( (/i(0)./i(r)} 
oscillates around zero for large r), while at low tempera- 
tures is a ferromagnet ((/z(0)./.t(r)) > for large r). How- 
ever, this scenario evolves at shorter distances through Tc 
where there is a strong change in the average character 
of the short-range magnetic interactions with the nearest 
neighbour (NN) and next near neighbours (NNN) ions. 
Focusing on the short ranged correlations in Fig. Ha it is 
apparent that, on approaching Tc from above for NN and 
NNN Mn-Mn pairs, FM correlations build up within the 
paramagnetic matrix. Notice, as shown in Fig. ^ that at 
T = 280 K small FM Bragg peaks start to develop, imply- 
ing the existence of already large FM domains (around 
100 A) before the system reaches the metallic regime. 
On lowering the temperature further, these FM domains 
grow to the point in which the basic background is com- 
pletely FM below Tc. Most interestingly, in addition 
to these FM correlations, we also observe that for those 
Mn-Mn pairs at the larger distances with respect to the 
Rietveld average distance (/i(0)./i(r)) gets into negative 
values around Tc (See Fig Hb). This suggests the ex- 
istence of magnetic correlations of an AFM character in 
this material with typical coherent lengths not larger that 
about 4-6 A. 

Figure |l| shows the normalized integrated intensities 
corresponding to the positive and negative lobes of the 
NN spin correlations since it is mostly for these that we 
have enough pairs (about 50 in the tails) to compute 
quantities that are sound on an statistical basis. Notice 
that the growth of a negative lobe in (^(0)./x(r)) does 
not correlate with the temperature dependence three di- 
mensional distortion of the octahedra derived from the 
Rietveld refinement, but it correlates well with the cal- 
culated from the Mn-0 pair distribution function of the 
RMC configurations (see Fig. Q). Furthermore, from the 
structural insight gained with the RMC, we can establish 
that, on average, the Mn pairs AFM coupled are located 
in the (a, c) plane (see Fig. |^), while those pairs at a 
distance closer to the average structure (Rietveld value) 
are along the long axis h. As a general trend, it is also ob- 
served that those Mn atoms exhibiting AFM interactions 
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are located at larger distances from the intermediate oxy- 
gen than its FM counterparts. 

In conclusion, our RMC analysis of the Bragg and dif- 
fuse neutron diffraction intensities draws an inhomoge- 
neous escenario in which three different phases coexist 
(FM and AFM domains within a PM matrix) in LCMO 
near Tmi- This microsegregation is probably driven by 
the distortion in the positions of the Mn atoms within 
the (a, c) plane that develops around Tc and Tmi- As 
a result, there are Mn atoms whose distance to the NN 
Mn atom is shorter than that corresponding to the av- 
eraged structure. These atoms present also shorter Mn- 
O distances and above Tc hold short range FM inter- 
actions within the PM matrix. As the temperature de- 
creases below Tc these FM clusters increase their coher- 
ence length and the matrix becomes FM. Also, within 
the (a, c) plane we have found a significant number of 
Mn atoms separated from their Mn NN a distance longer 
than the averaged Rietveld distance. These pairs (about 
15% of the NN at 280 K) interact antiferromagnetically 
around Tc and display longer Mn-0 distances making 
double exchange less efficient as overlap decreases. It is 
worth mentioning that the AFM domains show a phe- 
nomenology similar to that reported for the GE-state in 
Lao.sCoo.sMnOa near the Neel temperaturecll, making 
plausible to identify these islands as probable precursors 
of the CE-phasc. Notice also that the onset of the metal- 
lic state does not occur until the AFM correlations con- 
siderably weaken, pointing out that the AFM entities 
act as strong electron scatterers that add extra restric- 
tions to the electron transport in the PM phase. It is 
remarkable that in the corresponding Sr doped system, 
which remains metallic above Tc, local AFM cortejations 
are far weaker and {fi{0).fj,{r)) remains positiveL^I. It is 
realistic also to think that the application of an exter- 
nal field would inhibit the formation of this exotic AFM 
phase resulting in the decrease of the electron scatter- 
ing and the lowering the resistivity. The picture emerg- 
ing from our RMC analysis is, therefore, fully consistent 
with the most recent developments in the theory of man- 
ganites that postul^Sjequal density phase separation in- 
duced by disordeiauia. The simplified Ising model used 
by Moreo et al. contains the basic ingredients to ac- 
count for the observed PS around metal insulator tran- 
sitions predicted, but real systems seem to show an even 
richer phenomenology (three phases are experimentally 
observed versus the two phases predicted). Due to the 
small lengthscale of the observed FM and AFM domains, 
the question of whether these can be considered as proper 
thermodynamic phases is open to debate. In line with 
previous works, one could also invoke magnetic polarons 
as the distorting entities but then a reasonable hypoth- 
esis has to be added to the theory in order to account 
for the different short range correlations (FM as well as 
AF). Finally, we point out that our results also support 
the approaches based on more thermodynamic grounds 
that predict the presence of significant AF coujjlings and 
phase segregation for double exchange systemsEa. 
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FIG. 1. a): Temperature variation of the orthorombic dis- 
tortion in the RMC local and Rietveld variables, as defined 

in the text, b): Temperature dependence of the Magnetiza- 
tion (left axis) ;Normalized integrated intensities (right axis) of 
the FM lobe, A(-|-),(open circles) and of the AFM lobe, A(-), 
(solid circles) of the spin correlations of NN Mn-Mn pairs. 
Inset shows the first derivative of the measured resistivity of 
LCMO. 

FIG. 2. a) NPD pattern (circles) for T= 280K along with 
the corresponding residuals for the conventional Rietveld re- 
finement (middle profile) and the corresponding for our RMC 
approach (lower profile). b) Solid circles are the experimen- 
tal points, grey areas are RMC contribution for the magnetic 
scattering (diffuse + Bragg), dashed lines are for the nuclear 
diffuse scattering and black continuous lines are the sum of all 
RMC contributions( Magnetic diffuse. Magnetic Bragg, Nu- 
clear diffuse and Nuclear Bragg). c)A close up of the results 
at 260K pointing out the differences in the magnetic difltuse 
scattering corresponding to a non fitted configuration (thick 
line), obtained as explained in the text, along with the cor- 
responding contributions of the fitted configuration, symbols 
as in b). 

FIG. 3. Upper frame: Distribution function for Mn-Mn 
pairs, gMnMn{r). Lower frame: In (a, c)— plane (continu- 
ous line) and out of plane contributions (dashed line) to 
gMnMn-totai{r) for the Mn-Mn pairs. The partial distribu- 
tion fuction is now normalized to the Mn-Mn pair distri- 
bution function, computed without geometrical restrictions, 

gMnMn-total{r). 

FIG. 4. RMC spin correlation functions vs. distances be- 
tween Mn-Mn pairs, r, for all the temperatures explored. Up- 
per frame shows the long distance asymptotic limits, that 
coincide with the macroscopic magnetic behaviour. Lower 
frame: Detail of the spin correlations for NN and NNN man- 
ganese pairs. 
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